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We present an analysis of the electronic structure of perovskite-related iridates, 5d electron com¬ 
pounds where a subtle interplay between spin-orbit coupling, tetragonal distortions and electron 
correlations determines the electronic structure properties. We suggest via electronic structure cal¬ 
culations that a non-collinear calculation is required to obtain solutions close to the usually quoted 
ieff— 1/2 state to describe the t 2 g hole in the Ir^^id^ cation, while a collinear calculation yields a 
different solution, the hole is in a simpler xzjyz complex combination with a smaller LzjSz ratio. 
We describe what the implications of this are in terms of the electronic structure; surprisingly, both 
solutions barely differ in terms of their band structure, and are similar to the one obtained by a tight 
binding model involving orbitals with mean field interactions. We also analyze how the elec¬ 
tronic structure and magnetism evolve with strain, spin-orbit coupling strength and on-site Coulomb 
repulsion; we suggest the way the band structure gets modified and draw some comparisons with 
available experimental observations. 


I. INTRODUCTION 

A significant amount of work has been devel¬ 
oped recently on the interesting physical properties of 
perovskite-based iridates. These are Ir^+:d^ based com¬ 
pounds with structure related to that of the perovskite 
(metal cations surrounded by oxygen octahedra form 
the local environment) where dimensionality can be 
tuned by moving along the Ruddelsden-Popper series 
A^+iSr^Osn+i-^’^ Octahedral crystal field together with 
large spin-orbit coupling produce one hole in the t2^ 
band, forming aje//= 1/2 state that could somehow 
resemble the situation in cuprates: a single hole in a 
square lattice, with antiferromagnetic (AF) correlations 
dominating, but the t2^ manifold is involved in this 
case instead of the that occurs in cuprates. Based 
on these similarities, the appearance of superconductiv¬ 
ity has been speculated,^ and recently experimental ev¬ 
idences have started to emerge suggesting this could in¬ 
deed be the case,^’^ or at least that some Fermi sur¬ 
face features are very similar to underdoped cuprates.^ 
These systems are also interesting because the role of 
spin-orbit coupling in the electronic structure is not com¬ 
pletely clear, in particular it has been under debate what 
causes the insulating behavior in Sr2lr04, whether the 
system is a Mott or a Slater insulator.^ One of the 
keys to answer this question is to describe the electronic 
state the one t2^ hole is in. This has been described as 
a pure je//= 1/2 state^’^^ or as a mixture due to con¬ 
tributions from both je//= 3/2 and e^ states.To 
elucidate this issue, several iridates have been measured 
via different spectroscopic techniques, revealing that the 
L^/S^ ratio can be differentfrom the expected value 
of 4 that would occur for the pure je//= 1/2 state.At 
present it is not clear what the actual size of the different 
parameters involved in determining the electronic struc¬ 


ture should be (Hubbard U, spin-orbit coupling strength, 
band width, etc.), thus a careful study of the evolution 
of the electronic structure with these could shed light on 
the electronic structure properties of this family of com¬ 
pounds. 

Thin films of these iridates have been grown under dif¬ 
ferent conditions. With them, various experiments on 
the effects of strain in modifying their electronic struc¬ 
ture have been carried out. Both SrIrOs and Sr2lr04 
have been studied under various types of compressive and 
tensile strain. The results indicate that SrIrOs, being a 
paramagnetic metal in the bulk, can become insulating 
via both disorder-induced Anderson-type localization or 
by the application of strain.In Sr2lr04 a similar situa¬ 
tion occurs, the charge gap gets broadened by the appli¬ 
cation of compressive strain when compared to the bulk 
case,^^ but the X-ray absorption spectra (XAS) obtained 
in films grown on different substrates shows the opposite 
behavior,that the gap increases under tensile strain. 
There is a controversy here that we will try to address 
from electronic structure calculations. 

Sr2lr04 is a magnetic insulator that shows a pos¬ 
itive ©cw with a small ordered moment compared 
to the effective paramagnetic moment that has been 
measured.Some papers indicate the ordered moment 
can be produced via slight canting of the otherwise AF 
moments.^^”^^ The canting angle (see Fig. lb) occurs 
due to the deviation of the Ir-O-Ir angle from 180 ° via 
Dzyaloshinskii-Moriya interaction.It is to be under¬ 
stood if this canted AF solution is consistent with the 
positive ©cw observed.We will show ab initio from 
calculations including a non-collinear description of mag¬ 
netism that such a canting occurs naturally due to the 
structural distortion. 

We provide here an electronic structure analysis of 
both SrIrOs and Sr2lr04. We have analyzed the evolu- 
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FIG. 1. (Color online.) (a) Layered structure of Sr 2 lr 04 , 
where the oxygen octahedral environment around the Ir atoms 
is highlighted, (b) Sketch of an Ir-O-Ir bond, with the (ex¬ 
aggerated) canted magnetic order between Ir atoms, which 
arises due to the interplay of spin-orbit coupling and the non¬ 
linear hopping path through oxygens, (c) Schematic top view 
of the in-plane antiferromagnetic ordering of the Ir atoms. 
The net magnetic moments lie on the plane of the layered 
square lattice, slightly canted towards the neighboring oxy¬ 
gen, leading to a weak FM component. 


tion of their electronic structures with the application of 
strain, pressure, varying the on-site Coulomb repulsion, 
analyzing the role of spin-orbit coupling and comparing a 
non-collinear and a collinear description for magnetism. 
Surprisingly, we find that non-collinearity is crucial in 
determining the symmetry and character of the i 2 g hole, 
but the band structure of the system is barely affected by 
this orbital character. In the particular case of Sr 2 lr 04 , 
we show how a tight binding model captures many im¬ 
portant features of the behavior of the system, yielding 
also similar band structures even without including the 
canting of local moments or octahedral rotations. We 
have tried to elucidate in what limit of spin-orbit cou¬ 
pling vs band width and vs U the system is and how 
one can understand the properties of the ground state of 
these compounds. 


II. COMPUTATIONAL PROCEDURES 

Ab initio electronic structure calculations based on 
the density functional theory (DFT)^^’^^ have been 
performed using two all-electron full potential codes 
(wien2k^^ and Elk^^) on various Ir compounds, whose 
structure will be discussed below. The exchange- 
correlation term is parametrized depending on the case. 
We have used the generalized gradient approximation 
(GGA) in the Perdew-Burke-Ernzerhof^^ scheme for 
structural optimizations (atomic positions and volume 
optimizations). For treating on-site Coulomb repulsion 
in the 5d manifold, we have employed the local den¬ 
sity approximation+U (LDA+U) method in various ro¬ 
tational invariant flavors: in the so-called “fully local¬ 
ized limit” and for the elk calculations also with the 
Yukawa scheme.A non-collinear scheme for treating 
the magnetic moments was also used within the Elk 
code, and compared with the collinear formalism utilizing 
the same code. 

Regarding the wien2k calculations, calculations were 
performed with a converged k-mesh and a value of 
Ylrnt^max= 7.0. Spiu-orbit coupling (SOC) was intro¬ 
duced in a second variational manner using the scalar 
relativistic approximation.^^ The Rmt values used were 
in a.u.: 2.23 for Sr, 1.96 for Ir and 1.60 for O when study¬ 
ing Sr 2 lr 04 , and 2.50 for Sr, 2.00 for Ir and 1.63 for O 
when analyzing SrIrOs. 

We have carried out calculations in SrIrOs and 
Sr 2 lr 04 , using the structures from Refs. 39 and 40, re¬ 
spectively. We have simulated the effects of both tensile 
and compressive strains by fixing the a lattice parame¬ 
ter to that of several well-known systems typically used 
as substrates for thin film deposition (KTaOs: a= 3.989 
A, LSAT: a=3.868 A, SrTiOg: a= 3.905 A, LaSrA 104 : 
a= 3.755 A, MgO: a= 4.212 A) and relaxing both the c 
lattice parameter and the internal coordinates for each 
case. We have thus explored how the electronic band 
structure evolves under different degrees of strain and 
also analyzed the evolution of the electronic structure 
and magnetic properties as a function of U, the on-site 
Coulomb repulsion. 


III. IONIC PICTURE 

In all these iridates, we have Ir^+:d^ cations sitting in 
an octahedral environment with different degrees of dis¬ 
tortions. Because crystal-field splitting is larger than the 
Hand’s rule coupling strength, the ions are in a t^g con¬ 
figuration. Due to the large spin-orbit coupling typical 
in 5d electron electron systems like this, different eigen¬ 
states for the single t 2 ^ hole may occur. In the single-ion 
picture, in the presence of a tetragonal distortion that 
splits the t 2 ^ triplet, two situations may occur depend¬ 
ing on the relative strength of spin-orbit coupling with 
respect to the tetragonal distortion, namely: i) if spin- 
orbit coupling is just a perturbation, there will be an xy 
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singlet and the xzjyz doublet will be split by spin-orbit 
coupling into the xz ^ i yz (U eigenstates), or ii) if spin- 
orbit coupling dominates, one should see the splitting 
caused by spin-orbit coupling into a ]eff— 3/2 (four-fold 
degenerate) and 1/2 (two-fold degenerate, higher in en¬ 
ergy) states, that become further split when a tetragonal 
distortion is introduced as second order. The first situa¬ 
tion will lead to a spin-half ion with 1;^ = ± 1 {1jz/^z= 2), 
and the second situation would produce an Ir^+ cation 
in a state with the following expectation values: U = 2/3 
and S ;2 = 1/6 (L^ / S;^ =4). Thus, tracking the ratio of 
orbital to spin angular momenta, one can describe what 
particular ionic limit the cation is closer to. Typically, 
these systems are described as je// = 1/2 states,but 
in principle both situations (and any other intermedi¬ 
ate one) can occur, and they will, depending on different 
variables, such as epitaxial strain, pressure, etc. The 
question is what signatures this evolution may show in 
the band structure, magnetic properties or their strain 
dependence. 


IV. Sr2lr04 CALCULATIONS 

As we mentioned above, the analogy to the high-Tc 
cuprates and the recent hints appearing in the literature 
indicating possible superconductivity^’^ have drawn con¬ 
siderable attention to this system. The Ir^+:d^ cations 
sit in an elongated octahedral environment, where the 
t 2 g levels would be split by a local tetragonal distortion. 
This structure can be seen in Fig. la. The layered struc¬ 
ture leads to an in-plane square lattice of Ir atoms (Fig. 
Ic) with a singly unoccupied hole, somewhat similar to 
the situation in cuprates. In this case, the hole is in the 
t 2 g manifold (as opposed to the eg hole in cuprates) and 
we are dealing with more extended 5d electrons here, 
compared to the 3d in cuprates. We will see below in 
the band structure (Fig. 3) that the O p levels are away 
from the Fermi level in these iridates, which could pro¬ 
vide a difference with cuprates. This also occurs for the 
layered low-valence nickelates,^^'^^ which were also sug¬ 
gested as possible candidates for cuprate-like physics be¬ 
cause they show similar electron count on the Ni square 
lattice. AF coupling (with a slight canting) of the Ir 
moments is shown in Fig. Ic. Nearest-neighbor AF ex¬ 
change also resembles the in-plane checkerboard pattern 
typical in cuprates. 

In order to describe the electronic structure of the sys¬ 
tem, let us first try to draw some light on the controversy 
of how the band structure evolves with tensile strain. For 
this sake, we have performed ab initio calculations fixing 
the in-plane lattice parameter to various typical a values 
corresponding to the above mentioned usual substrates 
utilized to grow thin films on. For each of those, we have 
optimized the atomic positions and the c out-of-plane 
lattice parameter. 

In principle, one could reason that in the case of 
Sr 2 lr 04 , with a layered structure, increasing the in-plane 



FIG. 2. (Color online.) Evolution of the band gap of Sr 2 lr 04 
as a function of the in-plane lattice parameter (tensile [com¬ 
pressive] strain to the right [left]) for U= 3 eV in the “fully 
localized limit” utilizing a collinear description for the Ir mag¬ 
netic moments. We see that tensile strain leads to a larger gap 
by reducing the in-plane hopping. The trends are consistent 
for other U values but the actual metal-insulator transition is 
largely dependent on the U chosen. 


lattice parameter would tend to reduce the in-plane Ir-Ir 
hopping and the corresponding reduction in c brought 
about by tensile strain will not produce any additional 
band broadening. The latter makes sense because of the 
layered structure, as one can see in Fig. 1, where negli¬ 
gible direct Ir-Ir hopping along the c axis is anticipated. 
Thus, the expected situation from an electronic structure 
point of view would be an increased metallicity as it is 
compressed in the plane. This would be in agreement 
with XAS measurements in Ref. 22 and in disagreement 
with gap estimates as a function of strain in Ref. 21, that 
show the opposite trend. 

Figure 2 shows the results of our ab initio calculations 
for the evolution of the gap with strain for a particular U 
value (U= 3 eV). The trends shown would be consistent 
for other U values, but the particular gaps obtained, and 
also the point for the metal-insulator transition would 
differ (the evolution with U will be analyzed in further de¬ 
tail below in Section VII). We chose this (probably) large 
U value with a broad insulating region for the sake of il¬ 
lustrating more clearly the evolution with strain. We can 
see in Fig. 2 how the gap becomes larger as in-plane ten¬ 
sile strain is applied, as one would expect from the sim¬ 
plistic electronic structure arguments we just explained. 
As the material is compressed in the plane, larger band- 
widths are obtained and eventually the system becomes 
metallic at short enough Ir-Ir in-plane distances. 

We show in Fig. 3 the basic band structure of the sys¬ 
tem for two different substrates to illustrate further the 
way the gap gets smaller as in-plane compressive strain is 
applied. Again, calculations are presented for a large U 
of 3 eV so that a gap opens up and its strain dependence 
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FIG. 3. (Color online.) Evolution of the band structure of 
Sr 2 lr 04 with strain for U= 3 eV in the “fully localized limit” 
utilizing a collinear description for the Ir magnetic moments. 
Two band structures at different in-plane lattice parameters 
are presented: on the bottom a larger gap is found for 
3.989 A (as that of KTaOs) and on top, a smaller gap is 
obtained for a= 3.869 A (that of LSAT). We can see that the 
gap becomes reduced as further compressive strain is applied. 


can be analyzed more clearly. We see that the band clos¬ 
ing comes about without changing the main features of 
the band structure or the dispersions, it is very much a 
rigid band shift (plus some additional band broadening). 

We proceed to describe the electronic structure of 
Sr 2 lr 04 in more detail. As explained above, the Ir^+:d^ 
cation in a (distorted) octahedral environment will have 
one hole in the t 2 ^ manifold. Our ab initio calculations 
show that when the LDA+U method is used for any finite 
U, an in-plane AF ordering is obtained, and we have used 
this configuration for all the results analyzed throughout 
the paper. The relative magnitude of the spin and or¬ 
bital moments depends both on strain and the U value 
(the evolution with U will be analyzed in detail in Sec¬ 
tion VII). The effect of U is rather simple, increasing 
the on-site Coulomb repulsion leads to larger moments, 
both orbital and spin components. However, the effect 
of strain is not easy to predict a priori. 

We can observe in Fig. 4 the evolution of the Ijz/^z 
ratio as a function of strain for a fixed value of U= 3 eV, 
the same one we used for analyzing trends of the band gap 



FIG. 4. (Golor online.) Evolution of the L^/S^ ratio of 
Sr 2 lr 04 with strain for U= 3 eV in the “fully localized limit” 
utilizing a collinear description of the magnetic moments. We 
can see that the values are always close to 2.0, being closer 
to 2.0 for the more insulating solutions in the tensile strain 
limit. 


with strain. Let us recall that a value of this ratio closer 
to 2.0 corresponds to a hole in an ± i yz orbital, and 
values closer to 4.0 would in principle correspond to the 
so-called je//=l/2 solution. We observe that all values 
are closer to 2.0 except at the metallic limit, where the 
moments have almost vanished (a of LaSrA 104 ). In that 
case, the ratio is very large but it does not indicate a 
je//=l/2 solution was encountered. The tendency shows 
that the more insulating phases tend to stabilize a ratio 
closer to 2.0. 

Consensus in the literature describes almost unani¬ 
mously this material as being a realization ofaje//= 1/2, 
with a certain degree of mixing with the lower-lying je// = 
3/2 states. However, we have just seen that within a 
collinear formalism using wien2k (which is the method¬ 
ology utilized for all the calculations presented up to this 
point), no solution that resembles aje//= 1/2 state, in 
terms of Ijz/^z ratio (or the separate values of hz and 

expected for that simplistic single-ion solution) can be 
found. Changing U, the initial conditions for magnetism, 
or even the LDA+U flavor (around the mean field^^ was 
also tested) helps modifying slightly the actual values but 
not the main conclusions that the system is closer to hav¬ 
ing one hole in an xz ± i yz orbital. One could expect 
that for a spin-polarized system, the spin-up and down 
channels can be substantially separated in energy so that 
spin-orbit coupling can barely mix them. Yet, the purely 
ionic jeff =1/2 solution has an expectation value of Sz of 
only 1/6, maybe not enough to induce a large spin-up / 
spin-down Hand’s related splitting. However, this seems 
to be the case according to the calculations. 

It is worth noting again, as can be seen slightly in Fig. 
3, that changing the 1jz/^z ratio does not seem to mod¬ 
ify the band structure. In the cases presented there, the 
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FIG. 5. (Color online.) Projected density of states for Sr 2 lr 04 
in the limit of hz/^z— 2.0 (top, with a of MgAl 204 ) where the 
hole is clearly in an xzjyz combination with no xy character, 
and an increased ratio on the bottom (with a of LaSrA 104 ) 
where a larger xy character of the t 2 ^ hole is observed. These 
are collinear calculations with U= 3 eV within the “fully lo¬ 
calized limit”. 


application of strain basically leads to a rigid shift in 
the bands (plus some band broadenings due to increased 
in-plane hoppings), but the evolution to a larger ratio 
at more in-plane compressed lattices does not induce a 
change in the band dispersions. However, we will show 
now that a change in the band character of the i 2 g hole 
does occur even if the bands look unaffected. We will see 
that the 1jz/^z= 2.0 limit presents a hole that is fully 
composed of an xz ± i yz orbital, and going towards a 
larger ratio modifies the occupancies mixing different or¬ 
bital contributions. 

To analyze this change, we present in Fig. 5 two DOS 
curves at different in-plane strains. Even though we saw 
before that moving towards a larger in-plane strain barely 
changes the band structure (only a rigid shift), the evo¬ 
lution towards a more metallic state induces a larger xy 
mixing in the t 2 g hole. However, the je//=l/2 solution 
is not obtained, there is always much larger xzjyz con¬ 
tributions, which is incompatible with having a je//=l/2 
state. One could think that this has to do with the huge 
tetragonal distortion of the octahedral environment.^^ 
However, at large tensile strains, such as the case of set¬ 


ting a to that of MgAl 204 , the in-plane Ir-0 distance 
becomes equal (only 0.3 % difference) to the out-of-plane 
metal-anion bond length. Thus, the reason behind this 
unmixing does not have to do with the on-site energies of 
the different t 2 ^ levels that are split by the oxygen crystal 
field. Instead, it is the relative strength of spin-orbit cou¬ 
pling, that acts as a perturbation to the crystal field and 
splits the xzjyz orbitals forming 1;^ = ± 1 eigenstates,"^^ 
at least in a collinear calculation, since up to this point 
all the calculations presented were carried out within a 
collinear formalism to treat magnetism. It is worth not¬ 
ing that within this collinear scheme within wien2k, the 
Ijzj^z ratio is very much independent of the magneti¬ 
zation direction, whether it lies in the plane or out of 
the plane. The moments (both orbital and spin compo¬ 
nents) become larger (by up to 30%) when the magne¬ 
tization is out of the plane. We will describe in more 
detail the magnetization direction within a non-collinear 
description, which yields an orientation of the magnetic 
moments consistent with experimental observations. 

Calculations were also carried out with the elk code, 
both in the collinear and non-collinear formalism. Within 
both schemes, we take the initial magnetization to point 
in-plane along the Ir-Ir direction. As shown in Fig. 6b, 
the band structures within a collinear and a non-collinear 
scheme for magnetism are very similar, the same kind 
of dispersions and only some rigid shifts occur between 
both methodologies. When looking at the local magnetic 
order, we observe that the expectation values of L and S 
are not completely parallel, the ground state is a slightly 
canted antiferromagnet, produced due to the deviation of 
the Ir-O-Ir angle away from 180° (see the zoomed region 
in Fig. 6a). 

A small net moment is obtained in each Ir atom, apart 
from the AF component, yielding a net total angular mo¬ 
mentum (spin plus orbital) Jfe = 0.042. This arises due 
to the non-collinearity between the local magnetism of 
neighboring Ir atoms (see Fig. 6c,d). Within the non- 
collinear DFT calculations, the angle between the mag¬ 
netic moments in neighboring atoms is 151° (180° would 
be collinear antiferomagnetism). Moreover, the angle be¬ 
tween L and S is a = 8.1°. Therefore, the small net 
ferromagnetic (FM) component that arises due to the 
Ir-O-Ir angle, is expected to yield a weak FM signal in 
Sr 2 lr 04 . This is in qualitative agreement with exper¬ 
imental observations^^ showing an ordered moment of 
about 0.14 fiB- In our case, the value can be smaller be¬ 
cause only the value inside the muffin-tin spheres is con¬ 
sidered and the 5d electrons that cause it are very spread 
into the interstitial region. The ratios between orbital 
and spin components yield LjS = 5.4 and L^jSz = 5.6 
for this non-collinear scheme. These are even larger than 
the purely ionic description of the ]eff— 1/2 state, and 
are certainly very far from being a representation of a 
hole in an xz ± i yz state. However, this is not re¬ 
flected in the band structure, which barely changes when 
this ratio becomes so drastically modified. Remarkably, 
the band structures of the magnetic non-collinear solu- 
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FIG. 6. (Color online.) (a) Scheme of the Ir-O-Ir distor¬ 
tion in Sr 2 lr 04 , which gives rise to weak ferromagnetism, (b) 
Comparison of band structures calculated with a collinear and 
non-collinear scheme, (c) Vector plot in the Ir plane, showing 
the local magnetization for the different Ir atoms as obtained 
in the DFT non-collinear calculation, showing the canted an¬ 
tiferromagnetic state. Panel (d) shows a zoom close to the 
Ir atoms, showing a magnetization contour plot shaped as 
a t 2 g orbital as well as a small contribution from the oxy¬ 
gen atoms, (e) Sketch of a minimal model to describe the 
electronic structure of Sr 2 lr 04 , by taking into account near¬ 
est neighbors hoppings (only one t 2 g orbital is shown in each 
metal atom and the corresponding p^r orbital in the neighbor¬ 
ing oxygen). Red-blue orbitals are dxy orbitals, and green is 
an intermediate oxygen orbital, (f) Electronic band structure 
of the tight binding model, obtained by means of a Hartree 
Fock calculation of the Hubbard model. The parameters used 
are Xsoc — t and U — St. No Ir-O-Ir angle canting nor octa¬ 
hedral distortions have been taken into account in the model. 


tion are very similar to those obtained with a collinear 
scheme, as shown in Fig. 6b. For all these calculations, 
we have used a value of U= 2.7 eV, nevertheless in Sec¬ 
tion VII we will discuss the dependence of the ratio L^jSz 
with the on-site Coulomb repulsion U. 

Leaving apart the small canting of the Ir magnetic mo¬ 
ments, the main features of the electronic structure can 
be captured with a simple tight binding model for the t^g 


orbitals in a square lattice of the form 

H = Hnn + Hsoc + Hu (1) 

where TLatat is the nearest neighbors tight-binding hop¬ 
ping, Hsoc the SOC projected onto the t 2 g manifold 
and Hu the local electron-electron interaction which will 
give rise to magnetic order. Taking local t 2 g orbitals (la¬ 
beled by v^w) in each Ir atom (labeled by j), we only 
considered as non-vanishing hopping terms those involv¬ 
ing indirect hopping through the oxygen atom situated 
between them. This situation is sketched in Fig. 6c. 

HjSfN = (2) 

where Tyv,i,w,j can be easily obtained by symmetry con¬ 
siderations. 

We include SOC by calculating the matrix elements of 
the SOC operator for the t 2 g orbitals. 

Hsoc = ^socL • S (3) 

With such a nearest neighbor tight binding model, we 
have also introduced electronic interactions in a minimal 
way with an intraorbital Hubbard model 

Hu = Y^ (4) 

i^w 

and we solved the system with a Hartree Fock mean field 
approximation. 

With those ingredients, the following features are ob¬ 
tained. The ground state of the system at t^g filling is 
AF in agreement with the DFT calculations. Focus¬ 
ing on a particular Ir atom, the expectation value of 

the magnetic moment changes sign between the differ¬ 
ent t 2 g orbitals, but giving a net contribution, in agree¬ 
ment with the jeff = 1/2 picture. Let us recall that 
those jeff eigenstates look like the following in terms 
of the standard t 2 g orbitals: = 1/2; = ±1/2) = 

{\yz^ ±cr) T Act) T /v^? where a describes 

the spin. Magnetic order introduces the same kind of 
band splitting which ultimately gives rise to the in¬ 
sulating behavior. And finally, the band structures 
(Fig. 6f) resemble in dispersions and locations of max¬ 
imum/minima the ones obtained in our ab initio calcu¬ 
lations. Nevertheless, we have not implemented the Ir- 
O-Ir canting angle into this simple tight binding model, 
and therefore it is not able to capture the weak ferro¬ 
magnetism of Sr 2 lr 04 . Yet, the band structure obtained 
matches pretty well the main features of the electronic 
structure (see Fig. 6f vs b), and naturally produces a 
hole in a state close to the jeff = 1/2 ionic description. 


V. SrIrOs CALCULATIONS. 

We have also run calculations in the structurally sim¬ 
pler SrIrOs, whose low-temperature structure is an or¬ 
thorhombic distorted perovskite.^^ For imposing differ¬ 
ent in-plane strains simulating an epitaxial growth on 
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a lattice parameter (A) 

FIG. 7. (Color online.) Evolution with in-plane strain of the 
L^/S^ ratio for a non-collinear calculation in SrlrOa- Val¬ 
ues on the order of 4.0 and higher are obtained through a 
wide range of in-plane strains, particularly within the insulat¬ 
ing regime (whose extent depends on the particular U value). 
Results with U= 2.7 eV in the Yukawa-screened version are 
presented here, the trend with strain being systematic for a 
wide range of U values. 


a cubic substrate, we have used a tetragonal perovskite 
structure, allowing for octahedral rotations and relaxing 
these using the GGA-PBE exchange correlation poten¬ 
tial. Again, we have imposed different a lattice parame¬ 
ters as chosen above for the case of Sr 2 lr 04 , optimizing 
the atomic positions and the off-plane lattice constant. 

Experimental evidences from thin films grown on top 
of different substrates show that compressive strain ap¬ 
plied to SrIrOs increases the resistivity.^^ The situation 
in the perovskite-based system is different to the layered 
Sr 2 lr 04 . Here, compressing in the plane implies (due to 
the approximate volume conservation) elongating in the 
off-plane direction, so it is not clear how this can affect 
the overall band structure, and in particular the open¬ 
ing/closing of a band gap around the Fermi level. On 
the one hand, compressive strain increases the in-plane 
hoppings but it reduces the out-of-plane ones. The situ¬ 
ation is more complex here. 

To analyze this compound, we have carried out cal¬ 
culations using a fully relativistic non-collinear scheme 
using the elk code, setting up an AF coupling between 
Ir nearest neighbors (so-called G-type AF ordering) and 
introducing on-site Coulomb repulsion via the LDA+U 
method. By imposing different a lattice parameters and 
optimizing both the c lattice parameter and the inter¬ 
nal coordinates, we can obtain the evolution of the band 
gap with strain. Calculations show that the band gap 
becomes smaller as in-plane tensile strain is applied, in 
agreement with experiments. In particular, if a U= 2.7 
eV is used within the Yukawa-screened formalism (see 
Fig. 7), SrIrOs becomes metallic for a ^ 4.12 A, ap¬ 
proximately, but this crossing point depends heavily on 
the U value chosen. This trend is in principle consistent 



a lattice parameter (K) 

FIG. 8. (Golor online.) Evolution with in-plane strain of 
the L^/Sz ratio for a collinear calculation in SrIrOs- The 
evolution with strain is also towards a smaller ratio at larger 
tensile strain, but the values are in a completely different 
limit. These particular values were obtained for U= 3 eV in 
the “fully localized limit”. We see that, like in Sr 2 lr 04 , a 
solution with ratio close to 2.0 is the most stable one. 


with the experimental evidences described above. Elon¬ 
gating in the plane produces a larger compression along 
the c-axis in order to retain an optimal volume, and that 
small Ir-Ir off-plane distance eventually leads to a metal¬ 
lic bonding. We have not tested the situation with com¬ 
parable Ir-Ir distances in the plane, but we expect that 
for a sufficient compressive strain, the gap will start to 
decrease as well. 

Figure 7 shows the evolution of the Ijz/^z ratio ob¬ 
tained within a non-collinear scheme. We observe that 
the values are significantly higher than those we obtained 
within a collinear scheme using wien2k for Sr 2 lr 04 , and 
more consistent with those values yielded by Elk for 
Sr 2 lr 04 in a non-collinear calculation and with the hole 
description provided by the tight-binding model. These 
values resemble more the je//= 1/2 solution often quoted 
in literature to explain the electronic structure of this se¬ 
ries of materials, and in some cases the ratio even exceeds 
that of the single-ion solution. When the metallic region 
is reached, the ratio drops drastically. Again, the tran¬ 
sition point between insulating and metallic phase will 
depend strongly on the value of U chosen. We picked 
U= 2.7 eV (in the Yukawa-screened flavor) since it al¬ 
lows us to describe the trends more clearly due to the 
enhanced (probably unrealistic) gaps, and these trends 
with strain are consistent in a broad range of U values. 

Let us now try to draw some comparisons with our 
previous calculations in Sr 2 lr 04 using a collinear scheme 
as implemented in wien2k now for SrIrOs. The results 
are quite different, resembling more the picture described 
above for the layered compound. We can see them sum¬ 
marized in Fig. 8. We observe that again the values 
of the Ijz/^z ratio are closer to 2.0, just like we saw for 
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FIG. 9. (Color online.) Calculations for tetragonal SrIrOa 
with 1.5% uniaxial compressive deformation, (a) Comparison 
of the band structures with AF order (along the uniaxial axis) 
for both a collinear and a non-collinear scheme, (b) Evolu¬ 
tion with spin-orbit coupling strength of the L^/S^ ratio for a 
non-collinear and a collinear calculation. Evolution with SOC 
strength of Lz (c) and Sz (d), showing that switching on SOC 
triggers the magnetic order in the system. Noncollinear cal¬ 
culations yield larger Lz and Sz values, and also an increased 
LzjSz ratio. Calculations are performed with = 2.7 eV in 
the Yukawa-screened scheme. 


Sr 2 lr 04 . We even tried to introduce the density matrices 
of a solution as obtained in the non-collinear scheme in 
the collinear calculation, but it never converges without 
evolving towards the solution close to 2.0. This is quite 
unique, since the LDA+U method is usually capable of 
converging many different solutions, as long as the appro¬ 
priate density matrices are used at the initial step. The 
evolution with strain towards smaller ratios is consistent 
in both schemes, but on a completely different scale of 
values. 


from the non-relativistic limit \soc = 0 to its real value 
^soc = 1, where \soc is a multiplicative constant which 
scales the relativistic part of the Hamiltonian in the va¬ 
lence electrons 


H Hjion—relativistic T ^SOCH^p^iativistic (5) 

and we follow the evolution of the L^jSz ratio and the 
order parameters Lz and Sz- We perform such study for 
both a fully non-collinear scheme and a collinear scheme 
in the direction of the AF order. With this toy model, 
we obtained the Lz/Sz ratio as a function of spin-orbit 
coupling strength that we can see in Fig. 9 b. This shows 
an important difference in ratios only because of setting 
up a non-collinear calculation. As in the previous results, 
the collinear calculation yields smaller ratios throughout 
the different SOC strengths considered. 

Importantly, by following the evolution of Lz (Fig. 9c) 
and Sz (Fig. 9d), it is clearly observed that SOC trig¬ 
gers the magnetic order in the system. This can be un¬ 
derstood as follows: in the large SOC limit, the Fermi 
level lies in the middle of the jeff = 1/2 manifold, which 
has a smaller bandwidth than the full non-relativistic t 2 g 
manifold so that the system would have a large DOS at 
the Fermi level. In the intermediate regime, by switch¬ 
ing on SOC, the system undergoes a transition from the 
non-relativistic t 2 g limit to the relativistic jeff = 1/2 
state, which turns the systems vulnearable towards sym¬ 
metry breaking due to the decreasing bandwidths near 
the Fermi energy, in particular towards magnetism. 

Furthermore, non-collinear calculations yield impor¬ 
tant increases of Lz and Sz with respect to the collinear 
case. This difference also implies additional (minor) 
splittings in the band structure, as shown in Fig. 9a. 
It is worth to note that in the present system, the Ir- 
0-Ir bond is linear, so that this difference does not 
arise due to the canting of moments between different 
Ir atoms. Therefore, the difference in ratios between 
collinear and non-collinear schemes is related with the 
internal non-collinear nature of the magnetism within 
each atom. The present phenomenology suggests that 
non-collinearity might be of critical importance when ex¬ 
ploring jeff order. 


VI. EVOLUTION WITH SOC 

Trying to dig a bit more into the origin of this behavior, 
in this section we will study the dependence of the mag¬ 
netic/orbital order with SOC strength. We have set up 
a tetragonal SrIrOs structure with 1.5% tetragonal com¬ 
pressive uniaxial distortion in an AF order to be used 
as toy model. We stress that the incoming discussion 
is about a minimal model to understand the behavior 
of jeff within the collinear and non-collinear schemes. 
We performed LDA+U calculations with U = 2.7 eV in 
the Yukawa-screened scheme. In order to understand the 
effect of SOC, we have tuned the strength of the SOC 


VII. EVOLUTION WITH U 

In this section we will discuss the effect of the electron- 
electron interactions inside the d-manifold, as taken into 
account in an LDA+U scheme. In order to complete 
the previous discussion, we have performed fully non- 
collinear calculations with the elk package. For a sys¬ 
tematic study of the dependence with U of the electronic 
structure and magnetic properties, we have introduced U 
by means of the Yukawa scheme and fixing it to different 
values. The main results are summarized in Fig. 10. 

For Sr 2 lr 04 , magnetic order starts to show up from 
U = 1.5 eV, increasing rapidly with Lf ^ until reaching a 
value of S' = 0.06 for 7/ = 4 eV (which does not seem 
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FIG. 10. Evolution of the spin moment in the z direction 
for Sr 2 lr 04 (a) and SrIrOs (c) as a function of the on-site 
Coulomb repulsion. The panels on the right show the evo¬ 
lution of the ratio LzjSz for Sr 2 lr 04 (b) and SrIrOa (d). 
Both systems show an increasing ratio with increasing on¬ 
site Coulomb repulsion. All the calculations where performed 
within a non-collinear formalism. 


to be fully saturated). In comparison, the L^jSz ratio 
undergoes only small changes upon increasing U, showing 
a value around L^j = 5.7, above the expected value of 
4 for the ionic limit jf'e// = 1/2. 

For SrIrOs, however, (Figs. 10c and lOd) we can ob¬ 
serve a critical value of U from which magnetism sets in. 
It is worth noting that even when the system starts to 
become magnetic, it retains its metallicity until a larger 
U value is reached. For a lattice constant uq = 3.95 A 
the critical U for the onset of magnetism is Uc = 2.5 eV, 
whereas upon an artificial volumetric expansion of 5% in 
the lattice constant, the critical U lowers to Uc = 1-7 
eV. Even though such volumetric expansion is not fully 
realistic, it provides an insight on the effect of the Ir-Ir 
distance in the critical U. Regarding the ratio, it 


is observed that the ratio increases with U, but decreases 
upon increasing the lattice constant, as we explained in 
detail above when dealing with strain effects. 


VIII. CONCLUDING REMARKS. 

To summarize, in this paper we present ab initio cal¬ 
culations on Ir'^^-based oxides, which provide an elec¬ 
tronic structure understanding of the evolution of its elec¬ 
tronic structure, magnetic and transport properties with 
strain, some of them experimentally observed in detail. 
We also put on the discussion the issue of how accurate 
it is to describe these systems as je//= 1/2 states. To 
do this, we provide ab initio calculations based on differ¬ 
ent schemes together with an oversimplified tight-binding 
modeling of Srlr 204 . We can conclude that introducing 
non-collinearity effects on the Ir magnetic moments is a 
key ingredient in order to yield a solution that approaches 
the je//= 1/2 description of the t 2 ^ hole, as well as to 
capture the weak ferromagentic component. However, 
the effects on the actual band structure of introducing a 
solution that resembles the je//= 1/2 state are very mi¬ 
nor, and the evolution with strain is independent of how 
close to that solution the system is. A solution based on 
a combination of xzjyz orbitals plus spin-orbit coupling, 
which is obtained within a collinear description yields a 
very similar band structure and the same dependence on 
strain. However a non-collinear description is essential 
to determine the magnetic properties of the system and 
the orbital character of the single t 2 ^ unoccupied band. 
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